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We study the radiation of thermal photons and dilep- 
tons likely to be produced in relativistic heavy ion collisions. 
We find that the thermal photon multiplicity scales with the 
charged pion multiplicity as (dNch/dy)" with a ~ f.2 for a 
transversely expanding system, contrary to the general belief 
of a quadratic dependence. The scaling is shown to be valid, 
both for real and virtual photons. The coefficient of propor- 
tionality at a given energy may help us identify the appropri- 
ate equation of state of hot hadronic matter produced in such 
collisions. 

PACS numbers: 25.75.Dw, f2.38.Mh, 24.10. Nz, 25.75.Gz 



The search for the quark-gluon plasma (QGP), a 
deconfined state of strongly interacting matter, is en- 
tering a decisive phase. Evidence has been mounting 
for some time now, that such matter is perhaps be- 
ing produced in collisions involving heavy nuclei at the 
CERN Super Proton Synchrotron (SPS). If confirmed, 
this holds out a promise that the QGP likely to be cre- 
ated at the Brookhavcn Relativistic Heavy Ion Collider 
(RHIC) and the CERN Large Hadron Colhder (LHC) 
will occupy a much larger space-time volume, in similar 
collisions. The plasma thus produced will have to ex- 
pand and cool, due to large internal pressure. It has also 
to undergo hadronization, once the temperature drops 
below the transition temperature. These hadrons may 
continue to interact for some time before their density 
gets too low to support frequent collisions. We may also 
have a mixed phase of quarks, gluons, and hadrons dur- 
ing which the speed of sound may become vanishingly 
small, if quantum chromodynamics admits a first order 
phase transition. 

Many of the observables from the above collisions, 
like the particle distributions, electromagnetic radiation, 
the spatial and temporal dimensions of the fireball etc., 
will be decided by the equation of state of hadronic mat- 
ter. The transverse expansion of the system which dras- 
tically alters the time scales in the system, at least with 
regard to the mixed and the hadronic phases is also af- 
fected by the equation of state ||] . What should be the 
equation of state of hot hadronic matter which is created 
in heavy ion collisions? It is expected that if the tran- 



sition temperature is very close to the pion mass then 
the hadronic matter could possibly well be determined 
by pions and only a few low lying mesonic resonances. 
The hadronic matter is likely to be populated by heav- 
ier hadrons once the temperature is higher. Will these 
hadrons stay in chemical and thermal equilibrium, till 
the freeze-out? This will be decided by the competition 
between the collision frequency and the rate of expan- 
sion, and if valid, it will be reflected in uniquely deter- 
mined particle ratios g. It is not difficult to imagine 
that the equation of state or at least the composition 
of the hadronic matter may be somehow related to the 
manner in which the hadronization proceeds. This aspect 
has not yet received enough attention as it involves the 
non-perturbative realm of quantum chromodynamics. 

The question of hadronization of the plasma gets 
very complex in view of the recent findings that the QGP 
formed at RHIC or at LHC will not reach chemical equi- 
librium during the life-time of the quark phase [pi. We 
can understand this as follows. Recall that according to 
the Maxwell criterion, the phase transition temperature 
(Tc) is given by PgiTc) = Ph{Tc), where P stands for 
the pressure in the corresponding phase. The pressure 
for non-equilibrated QGP B can be approximated as 



PQ{T)^[mXg + 21X,]T^--B, 



(1) 



where A^ is the fugacity of the parton species k, B is the 
bag-pressure, and we have considered a two flavor plasma 
with zero baryochemical potential. Approximating the 
hadronic pressure in terms of a temperature dependent 
effective number of degrees of freedom Qh^T) [qh — 3, 
for a hadronic gas of mass-less pions, e.g.), we can write 



PH{T)=gH{T)-T\ 



which leads to 

rpi 



90B/7r2 
[I6A3 + 21A,] - gniT,) 
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We see that if we take the bag pressure to be a con- 
stant, then decreasing fugacities will, in general, increase 



the transition temperature. This is further accentuated 
by an increase in gu with temperature, which in fact re- 
flects the change in the composition of the hadronic mat- 
ter. The fugacities, measuring deviation from a chemi- 
cal equihbrium, may also have a radial dependence for a 
transversely expanding plasma |^ . These aspects encour- 
age us to think of a scenario where even the composition 
of the hadronic matter may become r -dependent. 

Till the development of such a complete description 
of the collision up to the point of freeze-out, hydrody- 
namics can provide a useful guideline. In as far as we 
can assume the validity of hydrodynamics, and addition- 
ally assume an adiabatic hadronization of the plasma, we 
can describe the collision in its entirety, provided that the 
equation of state of hadronic matter is known. 

In a recent study |j] we addressed the question of 
the equation of state in connection with the production of 
thermal photons and particles at CERN SPS in collisions 
involving lead nuclei, in a model calculation. Two ap- 
proximations were used to describe the hadronic medium. 
In the first one, only a small number of mesons (tt, p, w, 
and 77) was used ( " simplified equation of state" ) while 
in the second one we used all hadrons listed in the par- 
ticle data table g. Very different results were obtained 
for the transverse momentum distribution of photons and 
particles, depending on the assumptions made about the 
equation of state of hot hadronic matter as well as the na- 
ture of the initially created medium. In the present work, 
we show that a comparison of electromagnetic radiation 
for a number of final particle multiplicity densities may 
help us distinguish between the equations of state. We 
shall also see that the frozen-motion model, as a minimal 
extension |q| of Bjorken hydrodynamics H], provides us 
with a possibility to realize this distinction in an experi- 
ment done at a given energy. 

A clear illustration of the effect of the equation of 
state on the flow pattern at, say, LHC energies is seen 
in Fig. 1. We have plotted the space-time boundaries 
(at 2: = 0) of the QGP, the mixed and the hadronic 
phases, using the approach described in Ref. [^[||,^. We 
assume the system to be created in a state of QGP and 
use a boost invariant cylindrically symmetric transverse 
expansion. We see that the richer equation of state leads 
to a shorter-lived mixed phase and consequently the sys- 
tem edges past the freeze-out much more quickly. In fact 
this picture holds the key to the entire discussion that 
is to follow. One can, namely, immediately guess that 
the number of thermal photons or thermal dileptons pro- 
duced for the richer equation of state would be smaller. 
One may also witness a unique sensitivity on the extent 
of the strangeness equilibration |Iq] from such varied flow 
profiles of the system. These large differences could also 
be observed from pion interferometry. 

In order to proceed with our discussions we consider 
central collisions involving lead nuclei at a number of 
particle rapidity densities. Assuming an isentropic ex- 
pansion m of the plasma we relate the particle rapidity 
density idN/dy) to the likely initial temperature (Ti) as 



Tfn^ 



27r4 



dN 



45C(3) 4aQ7ri?2, dy 



(4) 



where aq = 377r^/90 for a QGP consisting of u and d 
quarks and gluons, Rt is the transverse dimension of the 
system, and we use the canonical value of 1 fm/c for the 
initial time r^. We shall fix the transition temperature at 
160 MeV and the freeze-out temperature at 100 MeV. 

The photon py-spectrum is obtained pT[ by convo- 
luting the rates for their production with the space-time 
history of the system as 



dN 
d^prdy 



- dr rdr d(j) drj 
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where the function fq (r, t) denotes the fraction of the 
quark gluon plasma in the system |l3] . The rates in the 
quark and the hadronic matter are evaluated according 
to methods developed earlier. Thus, we have included 
the Compton and the annihilation contributions from the 
quark matter, and the reactions involving tt, p, uj and 77 
mesons, along with the np —>■ ai —>■ ttj reaction in the 
hadronic matter |l3[ . 

We have shown the photon multiplicities as a func- 
tion of the multiplicity of the charged particles in Fig. 2. 
The slopes of the transverse momentum distributions for 
the two cases were found to be quite similar and are 
not shown here. We see that the "simplified equation 
of state" for the hadronic matter leads to about twice as 
many photons. This difference is likely to increase if the 
transition temperature is higher than 160 MeV, which we 
have used, as then the difference between the equations 
of state will become even more magnified. 

We have also marked the likely locations of the ex- 
pected values for colHsions at SPS, RHIC and LHC en- 
ergies using estimates for the particle multiplicities from 
Ref. |14| . It is quite clear that the availability of these 
results will go a long way in fixing the equation of state 
for the hadronic matter. From these results we can easily 
see that the thermal photon multiplicity is given by. 



dN^ 
dy 
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where a ~ 1.2 and K depends on the equation of state. 
In Fig. 3 we have repeated this analysis for the pro- 
duction of thermal dileptons, using the procedures dis- 
cussed in Ref. M , for dileptons having invariant mass M 
between 0.5 and 1 GeV. We have included the most dom- 
inant contribution of quark annihilation in the QGP and 
the pion annihilation proceeding via the formation of a p 
meson in the hadronic matter. We again find that twice 
as many thermal dileptons are produced for the simpli- 
fied equation of state (see also Ref. ^) and a scaling sim- 
ilar to Eq.(||) with a ~ 1.1 is satisfied. Apparently the 
slightly lower value of a in this case could be related to 



the fact that virtual photons carrying an invariant mass, 
are differently affected by the transverse flow. 

At first this slower increase of the number of thermal 
photons with the charged particle rapidity density looks 
very surprising, as it is generally believed that this depen- 
dence should be nearly quadratic for thermal sources [|5| . 
ft is not clear as to how this belief got currency, as it was 
pointed out by Feinberg ||l6[ two decades ago that the 
number of thermal photons or dileptons should scale as 
NJ , from very general arguments. In view of the im- 
portance attached to this scaling behavior for identifying 
new sources of photons and dileptons it is necessary to 
understand it more clearly. Consider a system consisting 
of A^ch charged particles. The number of thermal photons 
N-y will then be given by 



N^ 



e' iVch ly 



(7) 



where v is the number of collisions that each particle suf- 
fers. If the system lives long enough, as when it is con- 
fined to a box, every particle will have a chance to collide 
with every other particle, and then i/ '^ A^h- This will 
lead to the quadratic dependence suggested earlier p5[ . 
However, the number of collisions suffered by the parti- 
cles will be given by R/X, where R is the size and A is 
the mean free path of the particles, for a system created 
in heavy ion collisions |o] . Realizing that the number of 
particles will scale as W, we immediately get the scal- 
ing envisaged by Feinberg (see Ref. flfl] for arguments 
when the system may be undergoing expansion, which 
might bring in an additional factor of ■\/ln(Ach) ). We 
also realize that the transverse expansion of the system 
drastically reduces the time scales in the system 0. In 
the absence of the transverse expansion the life-time of 
the system can become extremely large; for example it 
will be several thousand fm/c, for the case discussed in 
Fig. 1 (see table I in Ref. 0). This would then mimic 
the case of particles contained in a box, and lead to the 
quadratic scaling reported by authors of Ref. ||l^] . 

Of course the transverse flow will considerably en- 
hance the production of photons having large transverse 
momenta, which however contribute only marginally to 
the total number. This is quite a sobering thought, as 
it suggests that it should be more difficult to disentangle 
thermal production from the background as the (total) 
yields in both cases show a very similar dependence on 
the final state multiplicity. The pt distribution of the 
background and the thermal photons will thus have to 
play a decisive role in identifying the latter. 

If one could assume that the minimal extension of 
the Bjorken hydrodynamics invoked by a number of au- 
thors PI to evaluate the production of thermal photons 
and dileptons at non-central rapidities is reliable, then 
one could also easily distinguish between different equa- 
tions of state by comparing the results at different rapidi- 
ties. It will be interesting to see whether a more detailed 
hydrodynamic model of the collision providing a fuller 
description at all rapidities |lq] can continue to support 



this distinction between the two equations of state. 

It is interesting to note that the excess dileptons 
observed at SPS energies for sulphur induced collisions 
by the CERES and HELIOS experiments |l|-||f|, which 
are measured at different rapidities, scale as (dNch/dy) , 
with f3 > 1. The statistics, however, is still not sufficient 
to distinguish the hadronic equation of state. 

In summary, we have evaluated the radiation of ther- 
mal photons and dileptons from relativistic heavy ion 
collisions for a range of multiplicities. Assuming the 
formation of a quark gluon plasma in the initial state, 
and properly accounting for the transverse expansion, 
we find that thermal production exhibits a much weaker 
than quadratic dependence on the final state multiplic- 
ity normally assumed. It is also shown that the results 
are quite sensitive to the equation of state used to de- 
scribe the hadronic matter. In order to distinguish and 
fix the equation of state we suggest to utilize the scaling 
rule of Eq. (^) , reasonably valid for both real and virtual 
photons, as the proportionality constant is determined 
by the equation of state. We are investigating whether 
this constant of proportionality is also a measure of the 
time-scale in the system Ref. |T^. 
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FIG. 2. Variation of rapidity density of thermal photons 
with the charge particle rapidity density for the two equations 
of state for the hadronic matter. 
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FIG. 3. Variation of rapidity density of thermal dileptons 
with the charge particle rapidity density for the two equations 
of state for the hadronic matter, having 0.5 < M 1 GeV. 
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FIG. 1. The space-time boundaries for the QGP phase, the 
mixed phase and the hadronic phase at LHC energies. The 
equation of state for the hadronic matter consists of only n, 
p, uj, and rj mesons (solid curves) and of full list of hadrons 
in the Particle Data Book (dashed curves). 



